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Abstract 

In this work, we explore the "degenerate gravitino" scenario where the mass difference between the 
gravitino and the hghtest MSSM particle is much smaller than the gravitino mass itself. In this 
case, the energy released in the decay of the next to lightest sypersymmetric particle (NLSP) is 
reduced. Consequently the cosmological and astrophysical constraints on the gravitino abundance, 
and hence on the reheating temperature, become softer than in the usual case. On the other hand, 
such small mass splittings generically imply a much longer lifetime for the NLSP. We find that, in 
the constrained MSSM (CMSSM), for neutralino LSP or NLSP, reheating temperatures compatible 
with thermal leptogenesis are reached for small splittings of order 10~^ GeV. While for stau NLSP, 
temperatures of Trh — 4 x 10^ GeV can be obtained even for splittings of order of tens of GeVs. 
This "degenerate gravitino" scenario offers a possible way out to the gravitino problem for thermal 
leptogenesis in supersymmetric theories. 

1 Introduction 

The existence of long-lived massive particles is always welcome in theories beyond the Standard 
Model (SM) if their lifetime is longer than the age of the universe (or completely stable) as they 
provide a candidate for the cold Dark Matter (DM) component of the universe to account for the 
abundance, QcDNih'^ — 0.11, inferred by cosmological observations. However, long-lived particles 
decaying after Big Bang Nucleosynthesis (BBN) place the theory in a very difficult situation when 
confronted with cosmological and astrophysical observation^. The very successful predictions of 
standard BBN are spoiled by the energetic products of the decay that can dissociate the produced 
light elements for lifetimes from 10^ to lO^'^ seconds. If the lifetime is between 10^*^ and 10^^ 



""^Decay products of particles with lifetimes as long as 10^ sec, see e.g. [ll[2j|3], can still produce observable 
signatures at present. Particles with longer lifetimes are only constrained from the measured relic density. 
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seconds, very stringent constraints come from the shape of the Cosmic Microwave Background 
(CMB) spectrum. For longer lifetimes, the emitted photons can reach us today and be observed 
as part of the Diffuse Extragalactic Background RAdiation (DEBRA). These observations restrict 
strongly the released energy in the decays of the long-lived particle at different times and hence its 
abundance, mass and lifetime. 

The paradigm of such long-lived particle conflicting with cosmological and astrophysical obser- 
vations is the gravitino in supersymmetric theories. The gravitino is the spin 3/2 supersymmetric 
partner of the graviton and the couplings of the gravitino to ordinary matter are gravitational cou- 
plings suppressed by the Planck mass, Mpi = (SttG^v)"^/^ ^ 2.4 X IQis GeV. Such small couplings 
make the typical gravitino lifetime (with electroweak scale masses and neglecting the masses of the 
decay products) of the order of 10^ seconds. Therefore, they decay after BBN and the gravitino 
abundance at BBN is severely constrained [1]. This has very important consequences in the phe- 
nomenology of this model. In fact, these BBN constraints forbid reheating temperatures larger 
than 10^ or 10^ GeV, which precludes thermal leptogenesis from generating a large enough baryon 
asymmetry. 

It is interesting to check whether it is possible to evade these stringent bounds. However, given 
that the gravitino couplings are completely fixed by supergravity, it is very difficult to change the 
gravitino production or its decay. The only "free" parameter available (from an effective theory 
point of view) is the gravitino mass itself and similarly the masses of the decay products. The 
gravitino gets a mass after supersymmetry breaking and the different soft-breaking terms receive a 
contribution proportional to the gravitino mass. If supergravity is the mechanism of mediation of 
the supersymmetry breaking from the hidden sector to the visible sector, we can expect the SUSY 
masses and the gravitino to be of the same ordei§|. However, most phenomenological analyses of 
BBN constraints assume that the released energy in the thermal plasma is of the same order as 
the gravitino mass itself. Although this is correct in most of the cases where the mass difference 
between the gravitino and the SUSY particles (typically the lightest MSSM supersymmetric particle 
(MLSP)) is sizeable, we can ask what happens if the gravitino and the MSSM LSP masses are much 
closer, AM = 7713/2 — "zmlsp ^ ?^3/2- This has two different consequences, first it is evident that 
the released energy in the thermal plasma is much smaller and this can help to relax the previous 
bounds. On the other hand, this small mass difference increases the gravitino lifetime and other 
constraints as the CMB spectrum or diffuse gamma rays come into play. 

In this work, we will analyze this scenario that we call "degenerate gravitino" scenario, including 
both the case where the gravitino is the NLSP decaying to the MSSM LSP and the case where the 
gravitino itself is the LSP and all SUSY particles decay into the gravitino. In the following, we 
define 6 as the degree of degeneracy between NLSP and LSP 

^ ^ "T-NLSP - "T-LSP _ AM 
"ILSP "ILSP 

In the degenerate gravitino scenario, the total dark matter abundance has two sources: one is 

^In principle, in gauge mediation or anomaly mediation models, constraints from gravitino decays can be 
evaded as gravitino mass can be, respectively, much smaller or much larger than typical SUSY masses 
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the usual LSP component from thermal production and decoupling and a second one from the NLSP 
non-thermal decays to LSP and SM particles. The sum of both components should reproduce the 
observed CDM relic density. The phenomenology of the degenerate gravitino scenario depends on 
both the identity of the lightest MSSM supersymmetric particle and on whether the gravitino is the 
LSP or the NLSP. The BBN, CMB and DEBRA bounds apply to the NLSP abundance while the 
LSP abundance is only constrained by the total dark matter abundance. In most cases, the lightest 
MSSM supersymmetric particle can be either the lightest neutralino or the lightest stau. Given that 
dark matter can not be a charged particle, we are left with three possibilities: gravitino LSP with 
neutralino or stau NLSP and gravitino NLSP with neutralino LSP. We will see that, for lifetimes 
smaller than the age of the universe, the strongest constraints on the reheating temperature arise in 
the case of gravitino NLSP, where all the previous constraints apply to the gravitino abundance and 
thus on the reheating temperature. On the other hand, when the gravitino is the LSP, it is possible 
to reach reheating temperatures compatible with thermal leptogenesis provided that the NLSP 
abundance at decoupling is sufficiently suppressed. Clearly, if the NLSP lifetime is much longer 
than the age of the universe, corresponding to small enough 6, the maximal reheating temperatures 
consistent with the observed DM abundance (Trh — 4.1 x 10^ GeV) can be reached. 

In this work, we will reanalyze the different constraints on the energy injected by NLSP decays 
and the reheating temperature in terms of this parameter 5 in the three above mentioned cases. In 
the next section, we will present the constraints on the energy release for different lifetimes of the 
NLSP in a model independent way. In Section [3] we apply these constraints to the gravitino case 
where the energy release and the lifetime are related. In section H] we analyse the case of the CMSSM 
looking for the new constraints on the reheating temperature in the degenerate gravitino scenario. 
Finally in section [5] we present our conclusions. 

2 Model-independent bounds 

If the NLSP lifetime is smaller than the age of the Universe, the energetic decay products can signifi- 
cantly affect cosmology. Even for lifetimes slightly larger than the age of the universe, observation of 
diffuse gamma rays can constraint this scenario. In this section, we consider the model-independent 
constraints on the degenerate NLSP-LSP scenario. As we will see, depending on the NLSP lifetime, 
these constraints originate from BBN observations, CMB spectral distortion or searches of diffuse 
gamma rays. There are additional constraints if the NLSP is charged. We focus on the case that 
the NLSP and LSP masses are nearly degenerate i.e. 5 <C 1 in Eq. ([1]). 

2.1 Relic abundance constraint 

The first constraint that should be considered is that the total relic density of cold dark matter 
must match the observed value. We will consider the situation where only the LSP and the NLSP 
are relevant; that is, heavier MSSM particles have already decayed to LSP/NLSP before BBN takes 
place. In general, the final abundance of the LSP will have two components. The first one is 
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the thermal abundance, ^^lsP' "^^ich comes from thermal processes occurring in the plasma like 
scatterings and freeze-out. The second one is the non-thermal component, ^^lsp*' that includes the 
contribution from LSP particles produced in NLSP decay. Therefore, the total cold dark matter 
relic density, will be % 

^CDM = ^^LSP + ^NLSP 1 (2) 

where the last term stands for f^LSp" ^'^'-^ have used Eq. ([1]). Assuming no late entropy release, 
Eq. ([2]) represents the present amount of dark matter, which should be equated with the observed 
value [S] 

^^WMAP = 0.1131 ± 0.0034 (3) 



It is useful to define the new parameter 



^ = T7 , (4) 

-iCDM 



which quantifies the amount of present cold dark matter coming from the NLSP decay. In this 
equation Inlsp refers to the NLSP yield just before its decay 0. Combining Eq. ([2]) and Eq. (jH) we 
get that 

^ = l_ "lsp ^ 

^WMAP h"^ ' 

which implies that u < 1 independently of 6. 



2.2 BBN constraints 

At temperatures of order T 1 MeV, the light nuclei are synthesized in the primordial plasma. These 
temperatures corresponds to times between 1 sec and 10'^ sec. The obtained abundances in standard 
BBN calculations are in striking agreement with observation 0. However, the injection of energetic 
particles in the primordial plasma at BBN or later can disrupt the standard BBN processes [51 El 
[TOl 111], leading to a disagreement between theory and observation. Thus, any additional particle 
decaying at BBN or later is subject to the strong constraints from light nuclei abundances. Usually, 
the energy of injected particles is assumed to be of the same order of the LSP/NLSP mass and 
stringent bounds were derived for various scenarios from BBN [T^fTHlfnifTOlfTUlfTTlfTHlfMl^mmi. 
These constraints are obtained by solving the full set of Boltzmann equations for BBN with a late- 
decaying particle. This requires the detailed study of the spectrum of decay products with all the 
relevant nuclear cross sections [9l \T0\ [TT] . The constraints apply to the released electromagnetic or 

■^Recall that we are assuming that the NLSP lifetime is smaller than the age of the Universe. If it were 
not the case, Eq. ^ is still valid as we are considering mNLSP — ™lsp- 

^Here the yield of a species i is defined as the ratio of the number density rii to entropy, Yi = rii/s. 
Recall also that the yield and the relic density of a massive particle species i arc related through Yi ~ 
4.1 X 10-12 (100 GeV/m,) (f], h'^/O.ll). 

^There are possible discrepancies in Lithium abundances, which might be explained using Gravitino dark 
matter [i[7]. 
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hadronic energy, parametrised by defined as [13] 



= B, Ynlsp , (6) 

where Ei is the released energy per decay with i = em for electromagnetic decays and i = had for 
the hadronic ones and Bi stands for the respective branching ratios. The constraints from hadronic 
processes are important when the lifetime is relatively short tnlsp ^ 10^ sec 116], The typical 
lifetimes considered in this work are larger that 10^ sec and therefore, in the following, we will 
consider only the BBN constraints on ^em- 

In the 2-body electromagnetic decay, the corresponding released energy is 

2 2 

TP "^NLSP ~ ^^LSP ('7\ 

-C'cm — 7, , \') 

^"iNLSP 

where we assumed that the visible particle mass is negligible. In the degenerate mass limit, 5^1, 
this reduces to E^.^ ~ ?71lsp ^- Notice that, as emphasized before, E'cm is much smaller than the 
usuall}!^ assumed value mNLSp/2. 

Given the above discussion, for our purposes, we can apply the constraints on ^cm from Ref. [11] 
directly to our scenario. Then, the constraint on reads 

Com ^ 4.1 X 10"i° Gev(^ ^^^^^Y^' ) ^ ^ < ^"PP- li-it ' 



where the right hand side in the above equation can be read off from the upper limit derived in [TT] 
for 5had = 0. 

For general values, the upper limit on the product oj i?cm ^ from BBN is shown in Figure [1] with 
red (dashed) line. From this figure, we can see that for lifetimes between 10^ and 10^*^ sec, this 
constraint requires LoBc^d to be between 10~^ and 10~^. Notice that BBN constraints apply equally 
to the case of charged NLSP decaying to gravitino and electromagnetic showers. 



2.3 CMB spectral distortion 

In addition to BBN constrains, for long lifetimes tnlsp <!/ 10^ sec, there are strong bounds from the 
shape of the CMB black-body spectrum. As pointed out in [12], the late injection of electromagnetic 
energy may distort the frequency dependence of the CMB spectrum from its observed blackbody 
shape. At late times of interest in our scenario, energetic photons from NLSP decays lose energy 
through processes such as — t- 7e~, but photon number remains conserved since other processes, 
like double Compton scatterings and thermal bremsstrahlung, become inefficient. As a result, the 
spectrum follows the Bose-Einstein distribution function 

^^^^^ ^ (.E/{kT)+^i _ I ' 

where /i here denotes the chemical potential. Then, the chemical potential of the distorted CMB 
spectrum has to satisfy the constraint, < 9 x 10~^ [23] . which, for decay lifetimes tnlsp ^ 

^Small mass differences has been considered in [3TJ[52] and for completely different motivations than ours. 
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Figure 1: Constraints on the combined parameters uBcm^ versus lifetime of NLSP using the result 
of [11] with Bh = 0. The red (dashed) line come from em BBN constraints and blue (solid) line from 
CMB distortion as explained in the text. Regions above the lines are excluded. 



8.8 X 10^ sec, translates into an upper bound on the released energy from NLSP decay, ^em defined 
in Eq. © [Ml [15], 

^em <1.59 X 10-8 e(-''c./mLSP)«/* ( 1^) GeV, (10) 

where Tdc — 6.085 x 10^ sec. 

For longer lifetimes (tnlsp ^ 8.8 x 10^ sec), the spectral distortions in the CMB spectrum can 
be described in terms of the integral of the fractional contributions to the energy e of the CMB per 
comoving volume during decay through the Compton y parameter, 4y = 6e/e, given by 

- = 7.04 X -^Cem, (11) 

where T{t) is the CMB temperature and t^s = [T{1 — /3)]^/^rNLSP, for a time-temperature relation 
T oc t~^, with F the usual Gamma function. In the radiation dominated era in the early Universe, 
for T < 0.1 MeV, 

/ t \-^/^ 

T = 1.15 X 10-3 GeV, (12) 

\ 1 sec / 

which gives /3 = 1/2. Thus teff = [F(1/2)]2tnlsp = vttnlsp- 

The observational limit \y\ < 1.2 x 10"^ [25] gives the constraint 



Com < 4.42 X 10-9 GeV (13) 

Tnlsp 
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Using Eqs. (|llll3p . this CMB constraint is plotted in Figure [T] with blue (solid) line. We can see 
that for lifetimes tnlsp ^ 10^*^ sec. the CMB constraint becomes more stringent than the BEN ones 
and sets the constraint on ^cm until recombination time. Similarly to the BBN, CMB constraints 
apply also to charged NLSP decaying to gravitino and electromagnetic showers. 



2.4 Diffuse Gamma-ray observation 

After recombination, at the cosmic time around 10^^ sec, the number density of free electrons drops 
quickly and the photons are almost free from the interactions. Therefore the photons from the 
decaying particles can reach us now and contribute to the Cosmic Gamma-ray Background (CGB). 
The observed CGB highly constrains any extra contribution including the photons from late decaying 
particles 0. 

The present photon flux from two-body decay can be written as 

c dt pc^WMApwS, 



dE^ Air 



1 dt p,OwMAPa;ijem _ (14) 

Ju ''"NLSP W-NLSP 



"iLSPi^ is the energy of the photon at production, tnlsp and mNLSP are the lifetime 
and mass of NLSP, pc = SH^/SttGn = 8.0992/i2 x 10"^^ GeV^ and a = a{t) is the time-dependent 
scale factor with a(to) = 1 at present time to- The delta function can be integrated using the 
formula 5{f{t)) = 6{t — ta)/\f'{t = ta)\ where ta is the solution which satisfies E^ = a{ta)Ecm- After 
integration, the flux reads 

d^ C Pcn^MAPCoBem e-W-NLSP 

dE^ 47r mNLSP TNLSP E^H{E^/Ecm) 

where the function simply cuts energies larger than the initial energy and each observed photon 
with energy E^ is produced at a time ta which satisfies a{ta) = E-y/E^^. Assuming the dark energy 
is a cosmological constant, this function ta{a) is given by 



t„.t(a = i.,/i..,J = 2Mv&^^M)ZVM, (10) 



where H{a) = Hq^/Qmcl ^ + ^^a- 

Taking into account that c/Hq = 1.3x10^^ cm, pc = 5.6x10"*^ GeV/cm^, Hq = 70 km sec~^ Mpc~\ 
Qa = 0.7, r^Af = ^CDM + = 0.3 and f^cDM = ^^wmap = 0.25, we find that the flux is 

"'^ 1.37 X 10^1 GeVcm-^ ^^^"^ ^ " =, (17) 



dE^ mNLSPTNLSP E^y/0.7 + 0.3{Ecm/E^, 



with 



ta = 3.51 X 10^^seclog[\/2.33a3 + y/l + 2.33a^]. (18) 



^Under appropriate conditions, the late decays of WIMPs to gravitinos and MeV photons may explain the 
MeV CGB anomalies 1221. 
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This differential flux must be compared to the observation of diffuse gamma ray flux. At each 
photon energy, we must require that this flux is smaller than the flux, E?^j^, observed by SPI, 
COMPTEL and EGRET [3]. Notice that, unlike the BBN and CMB constraints where we can find 
a bound on ^cm for a given value of tnlsPi now we need to specify both tnlsp and to obtain a 
bound on uBc^d. 

On the other hand, if the lifetime of NLSP is longer than the age of Universe, the line spectrum 
from the galactic center can be observed without cosmological redshift and this provides a further 
constraint on the emitted radiation from NLSP decays. For this we apply the bounds from Ref. [3]. 
Taking into account of the exponential decay of NLSP we have 



-e 



*°/"^^^^Cton<-^(^7 = ^cm), (19) 



47rmNLSP''"NLSP 

where psc = 0.3GeVcm~^ is the dark matter density at the solar distance from the Galactic center, 
Rsc = 8.5 kpc, Qim is a dimensionless integral of the line-of-sight intensity in the galactic center 
which ranges between 0.5 — 1.5 for various dark matter halo profiles [3]. The function is given in 
Figure 2 of Ref. [3]. Once again, as in the case of the CGB constraint in Eq. (jl7p . this constraint 
depends on the photon energy and therefore we need to specify both tnlsp and E.y to obtain bound 
on ujBpjri5. Thus, there is no simple analog of Figured] for a constraint on ^cm from Diffuse gamma 
rays observations. Moreover, the constraints on ^cm from the galactic center gamma rays are of 
the same order of magnitude (although slightly stronger) as the ones from the diffuse extragalactic 
emission. Therefore, for simplicity, we will only consider the diffuse gamma rays constraints, which 
apply to a broader range of energy. 



2.5 Catalyzed BBN 

Heavy long-lived negatively charged particles, X~, present during BBN can bind with light nuclei 
modifying standard BBN reactions. These catalyzed reactions, called CBBN, result in a change 
of light element abundances, and in particular lead to the overproduction of ^Li through CBBN 
reactions For lifetimes longer than 5 x 10^ sec, the observed light-nuclei abundances result on 

a constraint on the abundance of the charged relic, Yx~ < 2 x 10~^^ \27\ I28|. However, taking more 
conservative ^Li/^Li constraints, it is possible to relax slightly the previous bound to Yx- < 10~^^- 
10-15 [291 EQ]. 

Using Eq. (jU, we can translate the constraint on the yield of charged NLSP from catalyzed BBN 
to a bound on uj 

where we used Owmap^^ = 0.11 and Ycbbn is the maximum value allowed from catalyzed BBN for 
lifetimes larger than 10^ sec. As can be see from this equation, the catalyzed BBN bound is very 
stringent and indeed, it is very difficult to obtain such small yield at freeze-out in the MSSM |31j . 
However, it is still possible to find small allowed regions in the MSSM, or even in the CMSSM, where 
the f is the NLSP with large tan/3 [MIES] with relaxed ^Li/'^Li bounds [2T]. 
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3 Bounds on the degenerate gravitino scenario 



In this section, we will apply the model-independent bounds derived in the last section to the case 
of degenerate gravitino scenario. In our scenario, the gravitino can be either the LSP or the NLSP. 
In both cases it is thermally produced at reheating. After inflation, the Universe is reheated at a 
temperature Trh and gravitinos are produced through thermal scatterings in the plasma Their 
resulting relic density is linear in the reheating temperature and it is given by |35) 



OTP ,2 n 97 f A / 100GeV\ / A/3 V /01 ^ 

where M3 is the gluinos mass. Assuming no late entropy release, Eq. ^ sets an absolute bound on 
the reheating temperature, i.e. ilj^^ < Owmap, implies 

.„»0.v(^)(lf)\ (2.) 

This has to be compared with the minimum reheating temperature for successful thermal leptogenesis 
[37} [Ml [39] Trh ^ 2 X 10^ GeV. From this equation it is clear that the reheating temperature can 
not reach values much above ~ 10^" GeV. In addition to this constraint, we have to implement 
the constraints on the released energy considered in the last section. 

The relevant particles in the analysis of NLSP decays in the degenerate scenario are the gravitino 
and the lightest MSSM particle, which can be either the neutralino or the staiF^. We have two 
different situations depending on the particle nearly degenerate with the gravitino, namely gravitino- 
neutralino and gravitino-stau degeneracy. In each case the gravitino can be either the LSP or the 
NLSP. 



3.1 Gravitino-neutralino degeneracy 

As usual, in the MSSM, the lightest neutralino eigenstate Xi is parametrised as = Nii{—iB) + 
Ni2{—iW^) + Ni^H^ + NiiH^, where the unitary matrix defines the composition of neutralinos 
in terms of the Bino, Wino and Higgsinos. Since we are considering mass splittings that are smaller 
than the Z mass to suppress the hadronic branching ratio, the dominant (2-body) decay channel 
will be Xi ~^ 7 G or G — >• Xi7- As we will see below, the typical lifetime of NLSP is of order 
10^3 X (1 GeV/ AM) 3. Therefore, depending on AM we will have to consider different constraints. 

^We assume that at Trh, the Universe is composed of a thermal bath of MSSM degrees of freedom and 
that gravitinos production by inflaton decay is negligible jM] , 

^Taking into account the result of Ref. [36] the gravitino relic density would roughly increase a factor 2. 
possibility to increase Trh is to consider a squeezed gaugino spectrum reducing the gluino mass M3 
(see e.g. (40]). 

"'^"'^Sneutrinos LSP arc marginally allowed [41] due to their too large direct detection cross-sections. We do 
not consider them in our analysis. 
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3.1.1 Gravitino LSP 



In this case, the hfetime of neutrahnos is given by 

1.78 X 10^3 sec /iGeVV 



'^^^INiicosew + Nusmewl^y^Mj ■ ^^^^ 

Notice that, in the hmit of AM <C 7713/2, the hfetime depends only on the mass splitting 5 x 
7773/2 = AM, but not on the overall mass scale 7773/2 [H [l3]. If the gravitino is the LSP, the 
neutralino will decay into gravitino and photon and there are strong constraints on the released 
energy ^em = ^cm Different constraints will apply depending on the lifetime which, in turn is 
fixed by AM. 

• For 10 GeV < AM < 90 GeV, the relevant constraint is BBN and from Fig. [T]we have roughly 



d < 10-3. 



• For smaller splittings, from 1 GeV < AM < 10 GeV, the constraints from CMB spectrum are 
much stronger than the BBN ones and u6 is between 10"^ and 10"^ 

• For 30 MeV < AM < 1 GeV we have to take into account the diffuse gamma ray observations. 
Although the constrains on uj5 depend on 7773/2, for 7n,3/2 — 100 GeV, typical values range 
from uj6 = 10~^ to ujS = 10"^. 

• Mass differences from 2 MeV < AM < 30 MeV correspond to NLSP with lifetimes 2 x 
10^^ sec > tnlsp ^ 5 X 10^^ sec. These NLSP are already decaying at present, therefore 
diffuse gamma ray observations constrain their abundance. In this case uj6 ranges from 10~^ 
to 10"^ 

• Finally, for smaller AM, the neutralino is still present in the universe and it is completely 
stable for practical purposes. The only constraint comes from the WMAP measurement of 
the dark matter abundance. 

These constraints are summarized in Figure [2j We can see that for a given AM the constraints on 
ui are very strong for AM > 30 MeV, and, as we will see, it is difficult to reach such small neutralino 
abundances at decoupling in the MSSM. For smaller AM < 2 MeV, corresponding to neutralino 
lifetimes longer than the age of the universe, w ~ 1 is evidently allowed. For 2 MeV < AM < 30 MeV, 
uj ranges from 1 to 10"^. 

The only constraint on the gravitino abundance comes from Eq. ([2]) and the dark matter abun- 
dance measured by WMAP. Given u, we can use Eq. ([1]) and Eq. Q to calculate the required 
reheating temperature 

From Eq. (j24|) . we see that, provided uj <^ 1, one gets the maximal allowed reheating temperature 
Trh = 4.1 X 10^ GeV. 
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Q-6 llll I I I I I _ 

102 10" 10-' 10-2 10-3 10-4 

AM (GeV) 



'^3/2 — I Ar_ ^^^/i„ , Ar /!„ 12 ( A i\,r ) ' (^^) 



Figure 2: Constraints on the combined parameters W-Bcm versus AAd for given mass of neutralino 
NLSP (for A'^ii ~ 1) with Gravitino LSP. Sohd hne is for rriNLSP = 100 GeV and dashed hne for 
?7iNLSP = 200 GeV. The blue hnes for AM > 2 GeV come from BBN and CMB constraints and the 
green lines for AM < 2 GeV come from the diffuse gamma ray observations. Regions above the hnes 
are excluded. 

3.1.2 Gravitino NLSP 

If the gravitino is the NLSP, the dominant decay channel is G — t- 7 x5 and its hfetime is given by 

3.56 X 10^^ sec / I GeV " ^ 

lA^ii cos Bw + Ni2smew\'^ \ AM 

which is only a factor 2 larger than the gravitino LSP case. However, unhke the gravitino LSP case, 
the bounds on the released energy from this decay constrain strongly the initial thermal abundance 
of gravitinos. Notice that, here, oj represents the fraction of neutrahnos coming from gravitino decay 
and therefore the initial abundance of gravitinos. 

The constraints on uj5 for different ranges of AM seen in the previous section apply equally in 
this case. Now, given AM and 7/13/2, which fix 6, we have a direct constraint on the gravitino thermal 
abundance, and hence on u. As we know, the gravitino thermal abundance is directly proportional 
to Trh which in this case can be written 

Prom this equation, we see that, opposite to the case of gravitino LSP, in order to maximize the 
reheating temperature, one needs ui as large as possible. The constraints on Trh as a function of 5 
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2 Q4 UU I I I I I _ 

102 IQl 100 10-' 10-2 10-3 10-4 

AM (GeV) 

Figure 3: Constraints on the maximal reheating temperature after inflation for the Gravitino NLSP 
with neutrahno LSP (with Nn ~ 1) for rriNLSP = 100 GeV (sohd) and 200 GeV (dashed). Regions 
above the hnes are excluded. Notice that, as we can see from Eq. (j26|) . these maximal reheating 
temperatures are only possible if the corresponding bounds on cj are satisfied 

are shown in Figure [3l Given the strong constraints on u for the different 6s, Trh is considerably 
smaller than the maximal value allowed by WMAP. The only exception to this are the cases with very 
small AM ^ 2 x 10^'^ GeV, corresponding to the NLSP still present as a dark matter component, 
where w < 1 and the LSP abundance is smaller that the observed f^wMAP 

3.2 Gravitino-stau degeneracy 

In the MSSM, the lightest stau state ri can also be the LSP. However this situation is usually 
discarded as staus cannot play the role of cold dark matter. The remaining possibility is then that 
Ti is the NLSP with the gravitino as the LSP. The dominant decay channel in this case is ri — t- r G 
if the mass-difference is larger than the tau mass, i.e. AM > 1.77 GeV 0. In this case, the lifetime 

^^Notice that for AM < 1.77 GeV the two body flavour-conserving channel is closed and the stau can decay 
only through lepton flavour violating channels, ri — > /i G or ti — > e G, where the lifetime would be inversely 
proportional to the lepton-flavour violating coupling [321113]: t^:^ ~ 2 x 10^^ sec |^^,f^|~^ for m3/2 = 100 GeV 
and AM = 2 GeV. This means that, if these flavour violating couplings, the so-called Mass Insertions, are 
sizable, 6}^f^ > 0.03, the stau could decay before the present time. However, we do not consider this possibility 
in this paper. 
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of the stau NLSP is then given by 

.15 



r~ ~ 6 69 X 10- sec ( (^-^] ' ( ^-^ ^ (27) 



Notice that now, the stau hfetime, in contrast with the degenerate gravitino-neutrahno case, 
does depend on the overah mass scale ?n.3/2- Moreover for a mass difference similar to the tau mass, 
AM ~ 1.78 GeV, the stau lifetime is equal to the age of the universe. 

In principle, the stau has to satisfy similar constraints as in the degenerate gravitino-neutralino 
case with the exception of diffuse gamma ray constraints. Notice that the stau does not decay 
directly to photons and therefore does not contribute to the gamma ray background. Despite this 
fact, all the decays of the stau produce electromagnetic cascades that affect both BBN and CMB 
observables. 

Nevertheless, the main difference with the neutralino case is that stau, being charged, can form 
bound states with light elements and affect BBN predictions [281 [29]. Then, from Eq. ([20]) we obtain 
a strong constraint on u which for Ycbbn = 10~^^ is < 7.32 x 10~^mf /(300GeV). 

Finally, if the stau lifetime is longer than the age of the universe, the stau yield is very strongly 
bounded by the presence of anomalously heavy Hydrogen in deep sea water [H]. In terms of to the 
bound reads, w < 2.2 x 10^27 (m^/lOOGeV), for stau masses between 5 GeV and 1.7 TeV. Thus, in 
practice, this possibility can be completely discarded. In our analysis, we require that staus have 
already decayed at present and therefore we eliminate the staus with a lifetime longer than the age 
of the universe which corresponds to AM < 2 GeV. 

These constraints are presented in Figure [H where we see that for Ycbbn = 10~^^ the CBBN is 
the most stringent of all the constraints and requires a; < 7 x 10~^. This corresponds to stau relic 
density well below the observed relic density, Eq. ([3|), and the thermal relic density at freeze-out 
for purely RH staus [31]. In order to get such small abundance, the staus must have a substantial 
LR mixing. This requires both large /i and large tan/3 and moderate m^r^ [32^ I33| . which as we 
will see in section U is difficult, but still possible in the CMSSM. Furthermore, the vertical line at 
AM ~ 2 GeV corresponds to our requirement that all staus have already decayed at the present age. 

The bound on Trh here is analogous to the case of gravitino-neutralino degeneracy with gravitino 
LSP. Again, the reheating temperature is given in Eq. (j24p and given that, in this case, the allowed 
points in parameter space require lo < 10^^, one gets the maximal reheating temperature Trh = 
4.1 X 10^ GeV. 



4 The degenerate gravitino scenario in the CMSSM 

Perhaps the most appealing mechanism to transmit SUSY breaking from a hidden sector, where 
SUSY breaking occurs, to the visible sector is to use gravitational interactions, that are suppressed 
by the Planck scaled [l5]. This scenario is commonly called gravity mediation, and the CMSSM 

^■^Of course, in gravity mediated scenarios, there can be an associated cosmological moduli problem. In this 
work, as this issue is outside the focus of the paper, we do not address it and we will just assume that it is 
solved by some unspecified mechanism. 



13 




Figure 4: Constraints on uj versus AM with 77x3/2 = 300 GeV for stau NLSP with Gravitino LSP. 
The solid blue line represents the BBN and CMB constraints. The vertical line at AAf ~ 2 GeV 
corresponds to a stau lifetime equal to the age of the universe and therefore the region on the right 
is excluded from heavy-water searches. The dashed pink lines correspond to the CBBN constraint 
for ycBBN = 10-1^ and ycBBN = 10-1^ 

is one of its simplest and most popular realizations for phenomenological studies [l6]. It is defined 
in terms of only five free parameters: common scalar (jtiq), gaugino (777-1/2) ^-i^d tri-linear (^0) mass 
parameters (all specified at the GUT scale) plus the ratio of Higgs vacuum expectation values tan /3 
and sign(;u), where ^ is the Higgs/higgsino mass parameter whose square is computed from the 
conditions of radiative electroweak symmetry breaking. In addition to these parameters, in this 
analysis we have the gravitino mass, 7773/2, that we keep as a free parameter. 

The gravitino LSP (or NLSP) scenario in the CMSSM has been investigated thoroughly in the 
literature for the past years, both from the point of view of cosmological implications for dark 
matter [15] and for implications in collider searches [17]. One of the remarkable results of this 
scenario from cosmology is that the limits imposed on the reheating temperature of the Universe 
after an inflationary epoch have got down to a few 10^ GeVs [21]. This constraint basically rules 
out the thermal leptogenesis mechanism as the mean to produce the observed baryon asymmetry of 
the universe. 

In this section, we apply the results of the previous sections for the degenerate gravitino scenario 
to investigate the highest reheating temperature that can be reached in the CMSSM fulfilling the 
WMAP constraint on the CDM abundance and the relevant collider bounds; namely, direct SUSY 
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Observable 


Mean value 


Uncertainties 


ref. 






a (exper.) 


r (theor.) 




5a^ X 101° (e+e-) 


29.5 


8.8 


2.0 


m 


BR{B -> X,7) X 10^ 


3.52 


0.33 


0.3 


m 




0.1099 


0.0062 


O.ll^^/i^ 


m 




Limit (95% CL) 




r (theor.) 


ref. 


rrih 


> 114.4 GeV 




3 GeV 


m 


Sparticle masses 


As implemented in Micromegas 





Table 1: Summary of the observables used in the analysis to constrain the CMSSM parameter 
space. Upper part: Observables for which a positive measurement has been made. 5a^ = a^^^ — 
denotes the discrepancy between the experimental value and the SM prediction of the anomalous 
magnetic moment of the muon {g — 2)^. Lower part: Observables for which only limits currently 
exist. 

searches, {g — 2)^ using (e+e~ — T-hadrons) dat80 and the BR{B — > ^^7). All the constraints 
imposed are summarized in Table [TJ 

We use the fortran package SUSPECT [53] to solve the RGBs and to calculate the spectrum 
of physical sparticles and Higgs bosons, following the procedure outlined in [51]. To evaluate the 
CDM abundance in each point of the CMSSM parameter space, we employ the MicrOMEGAs code 
|52j . The branching ratio for the B — )• decay has been computed with the numerical code 
SusyBSG [53 using the full NLO QCD contributions, including the two-loop calculation of the gluino 
contributions presented in [SB] and the results of j57j for the remaining non-QCD tan /3-enhanced 
contributions. Finally we compute (^^^a^ at full one-loop level adding the logarithmic piece of the 
quantum electrodynamics two-loop calculation plus two-loop contributions from both stop-Higgs and 
chargino-stop/sbottom [58]. Then, the effective two- loop effect due to a shift in the muon Yukawa 
coupling proportional to tan^ /3 has been added as well [SS]. Recall that the communication among 
the different codes is done via the SLHA accord [60] . 

In particular, regarding the direct constraints on new particle searches, given that the theoretical 
error in computing the lightest Higgs mass ruh by SUSPECT is about 3 GeV [61], we require the 
calculated value of ruh to exceed 111 GeV. In the case of observables for which a positive measure- 
ment has been made, we require our predictions to be within the 2a range, for which we have added 
the theoretical and experimental errors, found in Tabled] in quadrature. In our numerical analysis 
we take nit = 173.1 GeV [62] . For the case of the stau NLSP, in addition, we completely exclude 
points that do not satisfy the conservative bound on 1~ < 10"^^ [29] from the catalyzed nuclear 
reactions [26] (see Fig. Hj) . 

We are now ready to present some representative numerical results. We are going to calculate 
the abundance of the MLSP (MSSM LSP), w = Vmlsp/VcdMi at freeze-out. Notice that w = a; if 

""^^Notice that we obtain the hadronic contribution to (g — 2)^ using only the data from (e+e^ — >hadrons), 
inclusion of the r data can decrease the discrepancy between the SM prediction and the experimental result. 
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Figure 5: Values of w in the ml/2-"^o plane of the CMSSM parameter space for nit = 173.1 GeV, 
^0 = 0, tan/3 = 10 and fi > 0. The red (medium grey) region is forbidden by the Higgs bound 
from LEP and in the very light grey region no correct electroweak symmetry breaking is obtained. 
The light brown (light grey) band below the NEWB area corresponds to the region forbidden by 
the LEP chargino bound. The regions below the dashed green lines satisfy the (g — 2)^ constraint 
at the 2 a (a^uSY > 2.43 x 10"^°) or 3 cj level (a^uSY > 1145 ^ IQ-IO). 



the MLSP is the NLSP, i.e. with gravitino LSP, but w = (1 — uj) when the gravitino is the NLSP. As 
a result, the reheating temperature will be proportional to (1 — w) in both cases. In the following 
we will focus on three representative tan/3 values: one moderate, tan/3 = 10, and two large values, 
tan /3 = 50 and tan /3 = 55, and we take Aq = 0. Indeed, we have performed several other scans 
varying both Aq and tan/3 and we found that there are no significant differences to these cases. 

4.1 Low-medium tan/3 = 10 

In Fig. [5]we show the region of w < 1 (corresponding to fi-^/i^ < 0.115) in the {mi/2, i^^o) plane for 
^0 = 0, tan /3 = 10, whereas we have chosen /i > motivated by (g — 2)^ data. In the CMSSM, 
the only mechanism which provides a dark matter abundance consistent with the WMAP value, 
if we require agreement with {g — 2)^ at 3 a, is stau-neutralino coannihilations |63j . This region 
is located in a narrow band above the stau-neutralino degeneracy line {rrir = rn^). In addition to 
this relic abundance constraint, the only effective constraint in this region comes from direct Higgs 
searches at LEP which excludes low gaugino masses mi/2 ^ 300 GeV. In this figure, the allowed 
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Figure 6: Maximal reheating temperatures in the mi/2-'mQ plane corresponding to the values of w 
in Fig. [5l 



values for w are shown in different colours. As we can see, in this case, we can only obtain w 
between 1 and 0(10"^). In the region of m^; < m^, although the annihilation mechanisms are more 
efficient than in the neutralino case due to fact that the stau is a charged particle, the lowest w 
we can reach is Wmin ~ 7 x 10"^, whereas from Eq. ()20p we see that the most conservative bound 
IV < 10~^^ from catalyzed nuclear reactions would require w < 2.7 X 10-2. Hence, this region of 

< rn^^ is completely excluded for all values of AM. If we eliminate the requirement of a non-zero 
SUSY contribution to {g — 2)^, we obtain two small regions. First, there is a vertical strip at high 
mo corresponding to a pole in the annihilation XX ~^ via the lightest Higgs. In this region we 
can obtain values of w as low as w ~ 0.01. Then, there is a long strip, below the non-electroweak 
symmetry breaking (NEWB) region, at very large mo, where the neutralino is a mixed Bino-Higgsino 
state and the annihilation to W'^W" is efficient enough to get a right relic abundance. However, in 
this region we obtain always w > 0.1. 

Using these allowed values for w and the corresponding u values, we obtain a bound on the 
mass difference AM in our degenerate gravitino scenario from Figure [2j According to Figure [21 for 
LO > 0.03 the relevant constraint is diffuse gamma ray observations, and this requires AM < 10"^ 
GeV. This means that the lifetime of the NLSP (neutralino or gravitino) is longer than the age of the 
universe. Notice that, in the case of gravitino LSP, the fact that we cannot obtain w = w < 0.1 implies 
that this scenario would be completely ruled out by BBN, CMB and diffuse gamma ray observations, 
unless AM < 10"^ GeV and the neutralinos are beginning to decay at present. However, if the 
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gravitino is the NLSP, we have two possibihties. The first possibihty would be that the observed 
dark matter abundance is provided completely by the neutralino w ~ 1, with only a small fraction, 
at the level of 10~^, due to gravitinos. This corresponds to Trh of the order of 10^ GeV, which is the 
usual situation in previous studies. The second possibility would be that there is a sizeable fraction 
of the dark matter due to gravitinos but, again, this would require Ai\/ < 10~^ GeV corresponding 
to a very long gravitino lifetime which permits to evade the cosmological bounds. The maximum 
values of Trh in these scenarios are shown in Fig. [6l As it can be seen from this figure, values of 
^RH > 10^ GeV consistent with thermal leptogenesis are accessible in the model, although this is 
only possible if AM < 10"^ GeV. 

4.2 Large tan f3 = 50 

In Fig. [71 we analyze the case of tan/3 = 50. For large values of tan/3, in the neutralino MSSM 
LSP region (m^^ < m^), we have different mechanisms to get a correct relic density consistent 
with WMAP. These mechanisms are i) stau-neutralino coannihilations, ii) "A-pole" region, where 
the s-channel exchange of the CP-odd Higgs boson. A, can become nearly resonant |64^ [65] and 
iii) the "focus point" or "hyperbolical branch" where a significant higgsino component, enhances its 
annihilation cross sections into final states containing gauge and/or Higgs bosons |66 [ I67 [ [5H] . In fact, 
the focus point region occurs at niQ much larger than mi/2 and therefore multi-TeV scalar masses 
which implies that 5a^^^^^ is well below the lower 2cr limit due to SUSY decoupling [69]. However, 
if we accept 

^a^issM ^Yid 3a level we can also reach the "focus point" region for mo ~ 2 x 10^ GeV. 
Both the neutralino-stau coannihilation band and the "A-pole" region merge for tan /3 = 50 as we 
can see in this figure. This region is cut from the small values by the lower 2a limit of the 

BR{B — )■ Xgj) constraint 0. This constraint limits the minimum possible value of w obtainable 
through the "A-pole". Therefore we can only reach Wmin ~ 0.1 both in the coannihilation and 
"A-pole" regions. In the "focus point" region the situation is again similar and we can reach only 
values Wmin ~ 0.1-0.01. 

In the case of rrir < m^, we can observe an allowed narrow band close to the tachyonic region 
in which the staus annihilate to bb via a pole in the s-channel exchanging a lightest Higgs. This 
process allows for values of Wmin ~ 10~^ which are consistent with the conservative CBBN constraint, 
1~ < 10~^'^-10~^^. The broader band at large mx/2 corresponds to the annihilation into a pair of 
light Higgses exchanging a light Higgs via a s-channel and here we obtain Wmin ~ 10~'^ and therefore 
1^: < 10~^^. On the other hand, we have to keep in mind that if we had imposed the more severe 
constraint 1^ < 10~^^, both regions would be ruled out. 

Likewise, with these values for w, we obtain, from Figure [21 the allowed values for the mass 
difference AM. In the neutralino-stau coannihilation or A-pole regions only AM < 10~^ GeV 
are allowed, similarly to the tan /3 = 10 case. However, in the region of < with values of 

^^The reason for this is the destructive interference of the chargino/squark loops which grows with tan/3 
[70] . This effect is enhanced for chargino masses < 100 GcV whereas is never dramatic for masses around 1 
TeV or larger. 
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Figure 7: Values of w in the mi/2-"^o plane of the CMSSM parameter space for mt = 173.1 GeV, 
^0 = 0, tan /3 = 50 and /i > 0. In this case, the red region represents the bound from BR(6 — t- S7) 
at 2a and there is a new region in dark brown (dark grey) marked "TACHYONS" corresponding to 
the presence of tachyonic masses. The other colored regions have the same meaning as in the case 
of tan (3 = 10. Regions below the dashed green lines satisfy the {g — 2)^ constraint at the 2 a or 3 
a level. 



w > 10~^, we can see from Fig. ^that values of the mass difference, 90 GeV > AM > 1 GeV, would 
be still allowed for these points lj. 

The maximum values of Trh are shown in Fig. [8l As we have seen before, both for the gravitino 
LSP and gravitino NLSP cases Trh (1 — w). The Maximal reheating temperatures for < 
are always Trh > 10^ GeV if the condition of AM < 10"^ GeV is satisfied, as in the case of 
tan /3 = 10. However, the region below the stau-neutralino degeneracy line are only allowed in the 
case of gravitino LSP when the maximal reheating temperatures are again Trh > 10^ GeV, although 
in this case we do not require a tight degeneracy between gravitino and stau. 

4.3 Large tan 13 = 55 

In Fig. [9l the case of tan /3 = 55 is explored. Similarly to the case of tan/3 = 50, in the neutralino 
MSSM LSP region, the main annihilation mechanisms are the neutralino-stau coannihilation and the 
"A-pole" resonance, however the "focus point" region is absent in this case. Again the BR{B — )■ X^')) 

^^We are considering only AAf < 90 GeV where the hadronic BBN constraints are not efficient. 
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Figure 8: Maximal reheating temperatures in the mi/2-"^o plane corresponding to the values of w 
in Fig. [71 Notice that the region below the stau-neutralino degeneracy line is not allowed in the case 
of gravitino NLSP. 



constraint limits the minimum possible value of w obtainable through the "A-pole". Therefore we 
can only reach w^-nin ~ 0.1, both in the coannihilation and "A-pole" regions. 

In the case of m;jr < ruy., as for tan /3 = 50, staus annihilate through an s-channel exchange of the 
light Higgs to a pair of light Higgses, although the annihilation to bb is suppressed in this case (it does 
not satisfy the required constraint Yf: < 10~^^). The allowed process, rr — )■ hh, allows for values of 
Wmin ~ 6 X 10~^ which are consistent with the conservative CBBN constraint for 1~ < 10~^^. 

Likewise, with these values of w, from Figure [2] we obtain the allowed values for the mass 
difference AM. As before, in the neutralino-stau coannihilation or A-pole regions only AM < 10~^ 
GeV are allowed and in the region of rUr < rriy. mass differences 90 GeV > AM > 2 GeV would be 
still allowed. The corresponding values of Trh are shown in Fig. [TDJ The discussion of section 14.21 
applies also in this case. 

5 Conclusions 

Even though they are attractive candidates for cold dark matter, gravitinos are usually a problem in 
standard cosmology. This is due to the fact that they typically decay at or just after the BBN putting 
in danger the successful light-elements abundances. Requiring that gravitinos do not conflict with 
big-bang nucleosynthesis implies that the reheating temperature should be low. In general, there 
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Figure 9: Values of w in the mi/2-'mo plane of the CMSSM parameter space for rrit = 173.1 GeV, 
^0 = 0) t^-n /3 = 55 and /x > 0. The meaning of the different regions is the same as in Fig. [71 The 
regions below the dashed green lines satisfy the {g — 2)^ constraint at the 2 o" or 3 o" level. 



is no experimental constraint on how low Trh should be, so in principle a reheating temperature 
as low as MeV, so to permit BBN, is perfectly allowed. However, successful thermal leptogenesis 
requires a high reheating temperature Trh ^ 2 x 10^ GeV. 

In this paper, we propose a solution that alleviates this tension by making the gravitino degen- 
erate with the lightest MSSM particle. This has the direct consequence that the injected energy is 
suppressed, making the decay products less dangerous for BBN. Due to the small mass splitting, the 
gravitino and the MSSM lightest particle are typically long-lived. We analysed this scenario (the 
"degenerate gravitino" scenario) by confronting it to cosmological and astrophysical constraints. 
Since the NLSP decays at or after BBN, we considered in addition to BBN, constraints from CMB 
spectral distortions and diffuse gamma rays observations. First we performed a model-independent 
analysis by considering a generic NLSP-LSP degenerate scenario where the NLSP decays through 
NLSP— )-LSP-|- X {X = 7,r). Since the final cold dark matter relic density is the sum of both 
thermal and non-thermal contributions, we required that the total cold dark matter is consistent 
with cosmological observation. Then, using the results of this analysis, we studied the degenerate 
gravitino scenario in the context of the CMSSM where three types of spectra arise, they are: grav- 
itino NLSP with neutralino LSP and gravitino LSP with neutralino or stau NLSP. Each of these 
cases has been analysed in this framework defined by the usual high energy parameters (mo, mi/2, 
Aq, tan/3, sign(^)) and the gravitino mass 7TI3/2, and confronted with low-energy observables. We 
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Figure 10: Maximal reheating temperatures in the mii2-mQ plane corresponding to the values of w 
in Fig. O As before the region below the stau-neutralino degeneracy line is not allowed in the case 
of gravitino NLSP. 



find that high reheating temperatures consistent with thermal leptogenesis can be found in all three 
scenarios if a sizable part of cold dark matter comes from gravitinos produced at reheating. In this 
case, depending on tan /3, we are led to regions in the parameter space where the relic density of 
stau and neutralinos are somewhat suppressed. In general to satisfy all the constraints, the mass 
splitting between the NLSP and LSP should be AM ~ 10~^ GeV for the degenerate neutralino- 
gravitino scenario, which implies very long-lived NLSPs which are beginning to decay at present. On 
the other hand, in the gravitino-stau scenario, splittings in the range 10 GeV < AM < 90 GeV are 
still consistent with reheating temperatures of the order of 10^ GeV if we consider the conservative 
CBBN constraint ycBBN < 10"^^ 

Let us comment on the required degeneracy in the "degenerate gravitino" scenario. Although a 
degeneracy of the order of AM ~ 10"^ GeV certainly implies a certain amount of fine-tuning, this 
tuning is only two orders of magnitude stronger than the usual tuning required in the coannihilation 
or funnel regions to obtain the right relic density in the MSSM. On the other hand, notice also that 
the fine tuning in our scenario is much softer that the tuning required in other scenarios like inelastic 
dark matter [7T] . 

Finally, it is also important to consider the phenomenological consequences of this scenario in 
colliders. In the case of neutralino LSP or NLSP, the only indirect signal of this scenario will be 
that the relic density of neutralinos inferred from the measurements of super symmetric masses and 
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couplings at LHC, will not match the observed cold dark matter abundance and will be smaller. How- 
ever, the measurements at direct detection experiments will agree with the cross sections obtained 
from colliders. On the other hand, for stau NLSP, the collider signatures would be spectacular, as 
the staus would be completely stable on collider scales and slow charged tracks will appear in the 
detector [72]. Notice that similar signatures can arise in other scenarios like for instance in dege- 
narate neutralino-stau scenario [l3] or in gauge-mediation scenarios. However, since the typical stau 
lifetimes in these scenarios are smaller than 1 sec, some of the staus will decay inside the detector. 
In contrast, stau lifetimes range from 10^ — 10^^ seconds in our scenario, making it very difficult to 
observe. Nevertheless, following the analysis of [73], stau lifetimes could be measured at LHC for 
mass splittings 30 GeV < AM < 90 GeV, corresponding to lifetimes 10 10 sec >Tr> 10^ sec. In this 
case, direct detection experiments will give a null results as all the dark matter at present times is 
made of gravitinos. Therefore, the "degenerate gravitino" scenario will be probed at colliders and 
direct detection experiments if SUSY is discovered at LHC. 

Note added: While completing this work we noticed the preprint [74j where the astrophysical 
consequences in dark matter halo properties of a neutral long lived decaying-dark matter particle were 
studied. The required parameters in our "degenerate gravitino" scenario in the case of neutralino 
LSP or NLSP are such that the constraints of [73] are satisfied. 
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